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Soft and semi-hard components structure in multiparticle production in 
high energy collisions 

A. Giovannini and R. Ugoccioni a * 

a Dipartimento di Fisica Teorica and I. N.F.N. - Sezione di Torino, 
via P. Giuria 1, 10125 Torino, Italy 

Superposition of jets of different flavors as well as superposition of jets of different topologies describe well observed 
structures (shoulder, H q oscillations) in e + e~ annihilation. The analysis of similar effects seen in experimental 
data in hadron-hadron collision is performed successfully by superimposing soft and semi-hard contributions at 
various energies. Negative Binomial multiplicity distributions are chosen in all examined classes of collisions as 
elementary substructures, i.e., as QCD inspired genuine self-similar fractal processes. Predictions on final particle 
multiplicity distributions in hadron-hadron collisions at 14 TeV are discussed. 



1. INTRODUCTION 

It is a fact, 1972 was a quite important year for 
multiparticle dynamics. It gave us 

1. QCD, the theory of strong interactions jl[, the 
non linear quantum field theory the past genera- 
tion was looking for in order to explain particle 
production firstly observed in cosmic ray physics 
in the thirties; 

2. KNO scaling behavior Q for final particle 
multiplicity distributions (MD's) which can be 
considered an extension of previous results at par- 
tonic level ||; 

3. the first evidence that final particle multi- 
plicity distributions in the accelerator region are 
broader than a Poisson distribution in clear dis- 
agreement with multi-peripheral model predic- 
tions [Q; 

4. the negative binomial (NB) distribution - 
called at the time Polya-Eggenberger || distri- 
bution — in order to take into account observed 
violations from multi-peripheralism in hadron- 
hadron collisions. 

How important all these discoveries were and 
still are in our field is witnessed by the fact that in 
our days, 25 years later, in approaching the study 
of the elementary substructures in multiparticle 



production in the new foreseen energy domain in 
hadron hadron reactions (see FELIX project at 
CERN LHC) all above mentioned results are to 
be used ||]. 

To disentangle elementary substructures in the 
new horizon is indeed the only hope we have in 
order to reach a simplified description of the ex- 
pected complexity of final states and to provide 
accordingly a better understanding of the subtle 
mechanisms controlling this extraordinary large 
number of final particles production and related 
correlations. 

Assuming our search successful, we have to 
answer at least two subsequent questions: are 
the above mentioned elementary substructures 
common to all classes of collisions? and — in 
addition — are they really elementary? 

Standard physical observables in this game are 
n charged particle multiplicity distributions, P n , 
factorial moments of the multiplicity distribu- 
tions, F q , and the corresponding cumulants mo- 
ments, K q , (these last observables are particu- 
larly sensitive to the tail of the distribution where 
events with many particles give a relevant contri- 
bution). Of particular interest turns out to be 
also the ratio K q /F q known in the literature as 
H q variable. 

All the mentioned observables are not indepen- 
dent, but they are linked by the following equa- 
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tions @,| 

OO 

F q = 5^»(n-l)"-(n-g + l)P„ (1) 

n— ij 

^=^-£( 9 7 V«-^ ( 2 ) 

Eqs. |l| and ^| — if interpreted correctly — show 
what is the main goal of multiparticlc produc- 
tion, i.e., the integrated description of n-particle 
multiplicity distributions and of the correspond- 
ing correlations. One should be able to relate the 
behavior of one observable to the germane behav- 
ior of the others; in addition the proposed expla- 
nation of one effect in one observable should also 
shed full light on the observed effect in the others. 

Now, both in e + e~ annihilation J?J and in 
hadron-hadron collisions (as we shall see) two ex- 
perimental facts on P n and H q observables in full 
phase space attracted our attention: 

a. the P n vs. n peculiar behavior (shoulder 
structure of n-particle multiplicity distributions) ; 

b. the H q vs. q oscillatory behavior (K q /F q 
g-particle correlations ratio oscillations). 

Here it will be shown that both effects can 
be understood in first approximation in hadron- 
hadron collisions in terms of the same cause, i.e., 
as the weighted superposition of soft and semi- 
hard contributions. This result should be con- 
fronted with the explanation of similar effects in 
e + e~ annihilation by the weighted superposition 
of two- and multi-jets contributions. 

Two remarks. Although the superimposed 
physical sub-structures in the two cases, e + e~ an- 
nihilation and hadron-hadron collisions, are dif- 
ferent the weighted superposition mechanism is 
the same. 

Secondly, all physical substructures are de- 
scribed by the same NB multiplicity distributions 
and corresponding correlation functions, which 
are QCD inspired genuine self-similar fractal pro- 
cesses. 



2. MULTIPLICITY DISTRIBUTIONS 
IN HADRON-HADRON COLLISIONS 
IN THE GeV REGION 

In the accelerator region it has been shown that 
final n charged particles multiplicity distribution 
in full phase space, P n , is initially narrower than 
a Poisson distribution then it becomes Poissonian 
and finally at higher cm. energies larger than a 
Poisson distribution. This behavior is usually de- 
scribed by using a two parameters distribution, 
the NBMD: the k parameter of the distribution 
is initially negative (the distribution is indeed a 
positive binomial) then it becomes infinite in cor- 
respondence of the the Poisson distribution and 
finally it is positive, i.e., it corresponds to a true 
NBMD. Now, being 1/fc closely related to the in- 
tegral over full phase space of the two particle cor- 
relation function above considerations are favor- 
ing anti-correlations in the lowest energy domain 
(k is negative, particles like to stay far apart), in- 
dependent particle production in the Poissonian 
regime and finally two particle correlation dom- 
inance (hierarchical correlations structure) when 
the multiplicity distribution is of NB type. 

The problem is that NB behavior for final 
charged particles multiplicity distributions can be 
trusted in hadron-hadron collisions in full phase 
space only up to ISR energies. At higher energies 
shoulder structures start to be clearly visible as 
shown by the UA5 Collaboration at CERN pp col- 
lider. The idea firstly suggested by C. Fuglesang 
is to explain observed NB regularity violations 
as the effect of the weighted superposition of soft 
events (events without mini-jets) and semi-hard 
events (events with mini-jets), the weight a so ft 
being the fraction of soft events and the multi- 
plicity distribution of each component being of 
NB type: 

■^n(^soft) ^-soft: ^soft j ^semi-harcb ^semi-hard) — 
a S oft-Pr B (™soft, fcsoft) + 

(1 - a sott )Pn B (n so{t , A: somi _hard) (3) 

In Fig. |l| are shown the results of the compo- 
sitions at 200 GeV, 560 GeV and 900 GeV cm. 
energies together with Pythia Monte Carlo cal- 
culations predictions; corresponding parameters 
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Figure 1. Charged particle multiplicity distribu- 
tions at cm. energies 200 GeV, 546 GeV and 900 
GeV. Data (white dots) are from UA5 Collabora- 
tion; the solid line is the fit with eq. j|, whose com- 
ponents are given by the dashed lines (the semi- 
hard component is shaded). The solid dots show 
the result of Monte Carlo calculations (Pythia) . 



Figure 2. H q vs q at cm. energies 200 GeV, 546 
GeV and 900 GeV. Data (dots) are from the MD's 
in Figure [l]; the solid line is the prediction of the 
fit also shown in Figure [l]. 
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Table 1 



Values of the parameters of eq. (||) fitted to UA5 
data fnj. 



c.m. energy 


200 GcV 


546 GcV 


900 GcV 


Qisoft 


0.93 


0.75 


0.72 


^soft 


19.9 


24.0 


26.9 


^soft 


7 


7 


7 


^semi-hard 


42.2 


47.6 


57.9 


^semi-hard 


79 


14 


13 



butions and correlations we emphasized in the in- 
troduction. What is really striking in our opinion 
is that elementary substructures of the weighted 
superposition mechanism in hadron - hadron col- 
lisions are again (as in e + e~ annihilation) NB 
multiplicity distributions. Analytic expressions of 
P n , F q and K q observables for the NBMD used 
in the text and in the figures are given in Table 

I 



Table 2 



NB observables used in the text. 



pNB 
r n 


k{k + l)...(k + n-l) 

n\ (n 


n n k k 
+ k) n+k 


pNB 


_ v (k + l)(k + 2)...(k + q 
= n H — : 

hi- 1 


-1) 


< B 


_,(<z-i)! 






(9-1)! 
(fc + !)(* + 2)... (fc + g- 


1) 



values of the two NB components are given in 
Table |l|. Notice that events with mini-jets (the 
shaded area in the figure) are highly depressed 
with respect to events without mini-jets at 200 
GeV c.m. energy: mini-jets contribution grows 
quickly going to higher c.m. energies. 

The conclusion is that Pythia Monte Carlo cal- 
culations predictions on P n are unsatisfactory in 
the GeV region, whereas the proposed fit in terms 
of the superposition of two NBMD's is quite good. 
It is interesting to remark that H q vs. q oscilla- 
tions (dots in Figure §) reconstructed from the 
n particle multiplicity distributions at different 
c.m. energies shown in Figure [l] do oscillate and 
that the first minimum is close to that observed 
in e + e~ annihilation. In addition by superimpos- 
ing weighted soft and semi-hard contributions H q 
vs. q oscillations are quite well described. 

In conclusion the same cause explains both 
shoulder effects and H q oscillations as requested 
by an integrated description of multiplicity distri- 



3. MULTIPLICITY DISTRIBUTIONS 
AND H q OSCILLATIONS IN 
HADRON-HADRON COLLISIONS IN 
THE TeV REGION 

Above mentioned results are important. The 
discovery that soft and semi-hard components can 
be described in terms of NBMD's allows indeed to 
model on purely phenomenological grounds and 
under very simple assumptions, expected scenar- 
ios for final particle multiplicity distributions and 
H q vs. q oscillations in hadron hadron collisions 
in the new energy domain to be opened at LHC. 
Accordingly we will study hadron-hadron colli- 
sions at 14 TeV c.m. energy, as expected by the 
LHC detector FELIX. 

The point is to find acceptable energy depen- 
dence of the NB parameters k and n for the two 
components substructures and the corresponding 
weight factor a so ft. In a region where QCD has 
no predictions we shall proceed by phenomeno- 
logical assumptions, which are discussed in the 
following. 

The first assumption concerns energy depen- 
dence of the total average charged particle mul- 
tiplicity, h. It is assumed here, as usually done 



n(y/a) = 3.01 - 0.4741nVi + 0.7541n 2 VJ (4) 

Since below 200 GeV c.m. energy one single NB 
fits multiplicity distribution data very well and 
above 200 GeV c.m. energy soft component has 
been disentangled as shown in Table |l|, we pro- 
pose to extrapolate the logarithmic increase with 
energy of the average charged particle multiplic- 
ity of the soft component, n so ft, also at higher 
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cm. energy according to the formula 

n soft (^ = -5.54 + 4.72 In y/i (5) 

For the average charged particle multiplicity of 
the semi-hard component, n se mi-hard> we then as- 
sume the UA1 result on mini-jets to be valid also 
at higher energy, i.e., 



^semi-hard (y/s) w 2n soft (v / s) 
It follows from 

fl = a so ftn S oft + (1 — OWt) ^semi-hard 

and Eqs. || and ||, that 



™ ^-semi-hard 

"soft = — = * 

^soft ^semi-hard 



n so {t 



(6) 



(7) 



(8) 



In this way the weight factor a S oft can be deter- 
mined at 14 TeV and it turns out to be equal 
to 0.30, indicating a large mini-jets and hard jets 
production. 

The second problem concerns the energy de- 
pendence of the second NB parameter, fc. In 
the accelerator region the parameter fc coincides 
with fcsoft; the soft component is indeed the dom- 
inant one. Above 200 GeV the global distri- 
bution is not of NB type and parameter k is 
hardly defined. However one can use the disper- 
sion D 2 = (n 2 ) — (n) 2 or equivalently the second 
factorial moment 



Fo D 2 - fi 
-4-1 = s— 



1 



(9) 



Indications from our fit for the soft component 
at 200 GeV, 560 GeV and 900 GeV cm. energy 
show that fc so ft is constant above 200 GeV cm. 
energy. In addition l/n so ft varies very little as 
the cm. energy increases. These two facts imply 
that KNO scaling should be valid for the soft com- 
ponent in a very wide energy range; we have no 
reasons to believe that this behavior will change 
in the TeV region. What about fc scm i-hard? Here 
two extreme scenarios are possible. 

Scenario 1 . KNO scaling holds also 

for events of the semi-hard component, i.e., 

-hard /^semi-hard = COnSt - This fact with the 



remark that l/n S emi-hard is a very small quantity 
in the TeV region leads to values of l/fc S emi-hard 
nearly constant but lower than l/fe so ft- Accord- 
ingly the resulting total fc value at 14 TeV is equal 
to 7.3 (fcscmi-hard is equal to 13, whereas fc S oft is 
equal to 7). The ratio D 2 /n 2 after a quick in- 
crease from the accelerator to ISR region reaches 
its maximum at approximately one TeV and then 
it decreases towards its KNO limit 0.16 at 14 TeV. 

Scenario 2. Strong KNO scaling violation 
for events of the semi-hard component, i.e., 



scmi-har 



d/ 71 : 



2 

scmi-har 



d is growing logarithmically as 



D 



D 

already shown for the semi-hard component in the 
cm. energy range 200 GeV - 900 GeV. We find 
fcscmi-hard to be equal to 3 and k total to 2.46 at 
14 TeV. 

Resulting n charged particles multiplicity dis- 
tributions at 14 TeV for the semi-hard and soft 
components together with the total multiplicity 
distribution for scenario 1 and 2 are shown in Fig- 
ures H and [| respectively. Notice that the mini- 
jets fraction of events is much larger than the soft 
events fraction. Expectations from Pythia Monte 
Carlo calculations at the same cm. energy are 
also given. The tail of the total multiplicity dis- 
tribution is highly suppressed in scenario 1 when 
compared to MD of scenario 2, suggesting that 
independent particle production is favored in the 
first with respect to the second scenario, where 
two particle correlations and larger mini-jets pro- 
duction are dominant. Pythia Monte Carlo cal- 
culations are closer to the soft component multi- 
plicity distribution in scenario 1 in the low mul- 
tiplicity region and to the semi-hard component 
multiplicity distribution in the large multiplicity 
region in scenario 2. Neither scenario reproduces 
Monte Carlo calculations for the total multiplic- 
ity distribution. 

In Figure [5] H q vs. q oscillations are also plot- 
ted in the two cases and fitted by H q values ob- 
tained via the weighed superposition of the two 
NBMD's. It is interesting to remark that in the 
KNO scaling violation limit oscillations disap- 
pear. 

One question should still be answered: are 
the substructures which we described in terms of 
NBMD's really elementary? The lesson which we 
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Figure 3. Charged particle multiplicity distributions at 14 TeV for scenario 1 (KNO scaling dominance). 
The long- and short- dashed lines represent the soft and the semi-hard component respectively (the latter 
is also shaded), the solid line represents the total MD. Dots are Pythia Monte Carlo predictions. 
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Figure 4. Same as Figure || but for scenario 2 (KNO scaling violation). 
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Figure 5. H q ratio vs q for the scenarios described in the text. The multiplicity distributions have not 
been truncated in order to calculate the moments. 



learned from residuals analysis of the 2-jet sample 
of events in e + e~ annihilation which was indeed 
quite well described by a single NB led us to the 
discovery that this sample of events was much 
better described by the weighted superposition of 
single jets of different flavor (each of them being 
understood as a NBMD), and to the prediction 
that single jets are dominated by two particle cor- 
relations, which are flavor independent. It seems 
therefore quite natural in order to answer to the 
above mentioned question to proceed firstly to 
the residual analysis of hadron-hadron multiplic- 
ity distributions in full phase space in the highest 
cm. energy region where data are available, i.e. 
again at 200 GeV, 560 GeV and 900 GeV cm. 
energies. 

Results of this analysis with one NB fit to the 
experimental data are shown in Figure |^a and 
with the weighted superposition of two NBMD's 
(each of them corresponding to the soft and to the 
semi-hard component respectively) in Figure ^b. 
The improvement from Figure ^|a to figure |^b is 
quite clear. This consideration notwithstanding 
we still see in the residuals of Figure ^b some ad- 



ditional substructures. When data on multiplic- 
ity distributions in full phase space will be avail- 
able at higher energies there are good reasons to 
believe that the situation will be very similar to 
that observed in the GeV region. 

Then the initial still unanswered question can 
be reformulated as follows: is it possible to asso- 
ciate the foreseen new elementary substructures 
in hadron-hadron collisions to single jets and to 
describe them — as done in e + e~ annihilation — 
in terms of NBMD's? 

Assuming the answer is yes, this fact would 
be really striking. It will teach us that a cru- 
cial step has been done in our understanding of 
NB universality. After this discovery NB univer- 
sality should be looked for not in the full sample 
of events, but at the most elementary level of in- 
vestigation — as suggested also by our studies on 
e + e~ annihilation — i.e., at single jet level in all 
classes of collisions . . . where one should find the 
domain of an effective self-similar Markov branch- 
ing process — as NBMD really is — at work. 

Assuming the answer is no, this fact would be 
even more interesting. It will tell us that we 
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Figure 6. Residuals for fits with (a) one single NBMD and (b) the weighted superposition of two NBMD's 
(shown in Figure @) at cm. energies of 200 GeV, 540 GeV and 900 GeV. 



reached in multiparticle dynamics a completely 
new domain whose fundamental mechanisms con- 
trolling particle production are all to be discov- 
ered. 



4. CONCLUSIONS 

H q vs q oscillations and P n vs n shoulder struc- 
ture in hadron-hadron collisions can be explained 
in terms of the same cause in the GeV region, i.e., 
the weighted superposition of elementary sub- 
structures which are identified with soft events 
(events without mini-jets) and semi-hard events 
(events with mini-jets) respectively. These results 
are then extended to the TeV cm. energy re- 
gion in two extreme scenarios characterized the 
first by KNO scaling behaviour and the second 
by KNO scaling violation of the semi-hard com- 
ponent. Extrapolations from both scenarios dis- 
agree with Pythia Monte Carlo calculation pre- 
dictions. It is interesting to remark that H q vs. 
q oscillations disappear in the KNO scaling vio- 
lation framework. A large mini-jets and jets pro- 
duction is expected in both scenarios. 70 percent 
of events are indeed of this type: this number 



is one order of magnitude larger than the cor- 
responding number at 200 GeV cm. energy. In 
addition produced mini-jets in the KNO scaling 
violation framework are expected to contain more 
particles than those in the KNO scaling scenario. 

These results are particularly interesting in 
view of the following remarks. 

1. The same weighted superposition mechanism 
of two components explains both effects also in 
e + e~ annihilation. 

2. Above mentioned substructures (soft and 
semi-hard events contributions in hadron-hadron 
collisions and 2- and multi-jets events contribu- 
tions in e + e~ annihilation) are all well described 
in terms of NBMD's. 

3. Residuals analysis of the 2-jet sample of events 
in e + e~ annihilation reveals further substructures 
which we associate to single jets of different flavor 
and describe again in terms of NBMD's. 

4. Residuals analysis in hadron-hadron collisions 
multiplicity distributions extrapolated in the TeV 
region reveals further substructures also in the 
soft and semi-hard components, whose study is 
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a challenging problem for future hadron hadron 
detectors. 

The next step of our programme is to extend 
the present theoretical investigation to rapidity 
intervals and impact parameter space. 
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